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Abstract 

Diffraction is one of the universal phenomena of physics, and a way to overcome it has always 
represented a challenge for physicists. In order to control diffraction, the study of structured waves 
has become decisive. Here, we present nondiffracting spatially accelerating solutions of the Maxwell 
equations: the Weber waves. These nonparaxial waves propagate along a parabolic trajectory 
while preserving its shape to a good approximation. They are expressed in analytic closed form 
and naturally separate in forward and backward propagation. We show that the Weber waves are 
self-healing, can form periodic breather waves, and have a well-defined conserved quantity: the 
parabolic momentum. We find that our Weber waves for moderate to large values of the parabolic 
momenta can be described by a modulated Airy function. Because the Weber waves are exact 
time-harmonic solution of the wave equation, they have implications to many linear wave systems 
in nature, ranging from acoustic and elastic waves to surface waves in fluids and membranes. 

PACS numbers: 42.25.-p, 03.50.De, 41.20.Jb, 41.85.-p 



An accelerating beam is a shape-preserving light beam whose peak intensity follows a 
continuous parabolic curve as it propagates in free space. The first accelerating beam, the 
Airy beam, was proposed and observed in 2007 pQ. By now, it has been proved that the 
only solutions to the paraxial wave equation with accelerating nondiflractive properties are 
the one-dimensional Airy beams and the two-dimensional Airy and accelerating parabolic 
beams [21 [3] . The unique properties of the accelerating beams has led to many intriguing 
ideas and applications; e.g., particles and cells micromanipulation [HE], plasma physics [6], 
nonlinear optics [7], plasmonics [HIE], and micromachining |10j . among other. 

Recently Kaminer et al. in [TT] introduced new nondiflractive accelerating waves called 
"half a Bessel" waves. These waves propagate along a circular trajectory, during a quarter 
of the circle they are quasi shape-preserving and after this diffraction broadening takes 
over and the waves spread out. The importance of these waves consist in having the same 
characteristics as the paraxial accelerating beams [IH3] but in the nonparaxial regime, i.e., 
these waves can bend to broader angles. Therefore, the "half a Bessel" waves allow to 
extrapolate all the intriguing applications of accelerating beams to the nonparaxial regime, 
and because these waves are solutions to the wave equation, they have implications to many 
linear wave systems in nature, ranging from acoustic and elastic waves to surface waves 
in fluids and membranes. Therefore, a natural question is, are there other accelerating 
nondiflractive solutions to the Maxwell equations? 

In this Letter, we present new nonparaxial spatially accelerating shape-preserving waves: 
the Weber waves. These nonparaxial waves propagate along a parabolic trajectory while pre- 
serving its shape to a good approximation during long distances (many diffraction lengths). 
Our Weber waves, as the "half a Bessel" waves, are self-healing, they can form breather 
waves, and they are a complete and orthogonal family of waves. However, the Weber waves 
naturally separate in forward and backward propagation, and they have an analytic closed 
form solution. 

The "half a Bessel" waves by construction break the circular symmetry and do not have 
well defined angular momentum. In contrast, we show that the Weber waves have a well 
defined conserved quantity: the parabolic momentum. We found that our Weber waves for 
moderate to large values of the parabolic momenta can be described by a modulated Airy 
function. We demonstrated that the local value of the linear momentum density of the 
Weber waves flows in parabolic lines and this allows to transfer parabolic momentum to 
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mechanical systems. Accordingly, schemes and uses equivalent to those found for orbital 
angular momentum [12J can be proposed for parabolic momentum. 

To construct new nondiffr acting accelerating waves of the Maxwell's equations as in 
one just need to find quasi-nondiffracting accelerating solutions of the two dimensional 
Helmholtz equation (d xx + d zz + k 2 ) ip — 0. Because the waves of interest are nonsingular 
and we consider only propagation in free space, there can be no evanescent components. 
Therefore, the spectrum of all the solutions to this equation lies on a circle or radius k = u/c 
in the k x — k y plane and it is characterized by the spectral function A (9) = {A + (9) , A_ (9)}. 
From this spectral function, the wave can be constructed by 

ik(x cos 9-\-z sin 8)^Q ^-Q 

We want forward propagating waves in the z-coordinate, therefore we are interested only 
in solutions with components in the top half < 9 < it of the spectral function, i.e., 
v4_ (9) = 0. In the rest of this Letter we will assume forward propagating waves; backward 
propagating waves can be readily constructed by A (9) — > {0, A + (—9)}. By this argument, 
the "half a Bessel" waves are defined in [11] as half of the spectral function of a Bessel beam, 
i.e., A + (9) = exp(im6) / y/n . The "half a Bessel" waves do not have analytic solution in the 
x — z space because Eq. Q for its spectral function cannot be solved analytically. 

Because we want solutions that have a parabolic translation upon propagation, we look 
at the two-dimensional Helmholtz equation in parabolic coordinates x + iz = (rj + i^) 2 /2 
with ranges in rj G (— oo, oo), £ £ [0, oo). In [13] a solution whose spectral function just lies 
on the top half of the spectral circle is created by superposing two parabolic nondiffractive 
beams. There, this solution is considered as a nondiffractive wave that propagates in a 
straight line, here we will consider that the propagation direction is transverse to our beam 
profile and not perpendicular. As it is show in [TT] this simply physical reinterpretation of 
the solution has profound implications. Therefore we define the Weber waves as [13], [H] 

|ri| 2 P (o-£;a)Po(cT?7;-a) ^ 

+2*|r 3 | 2 PiK;a)P 1 (cT77;-a) J ' 

where Ti = T(l/4 + ia/2), T3 = T(3/4 + ia/2), a = y2k, a is a dimensionless parameter of 
the solution that we will call the parabolic momentum and P 7 (•) are the parabolic cylinder 
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FIG. 1. (Color online) Top row (a-c) and bottom row (g-i): "Half a Bessel" waves; middle row 
(d-f): Weber waves. The dashed black parabolas describe the trajectory of the maximum lobe of 
the Weber wave in the same column. All figures depict a 220£;~ 1 x 220£;~ 1 square. 



functions [15] that in terms of the Kummer confluent hypergeometric function iF\ 
be writen as 



P 7 (x, a) 
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The spectrum of the Weber waves can be obtained from 
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For two dimensional waves the physical inner product, i.e., invariant under Euclidean 
transformations, is given by (VhlV^) — fl A\ (9) A 2 (9) d9 [16J. With this inner product, 
the Weber waves are orthogonal with respect to the parabolic momentum, i.e., 



{W( V} Z } a)\W( V} ^a')}=6(a-a>) 
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where 5 (•) is the Dirac delta function. Therefore, the Weber waves form a complete and 
orthogonal family of forward propagating waves and any forward propagating wave can be 
represented as a linear superposition of Weber waves with different parabolic momentum. 

We will now analyze the physical characteristics of our Weber waves. The middle row 
of Fig. Q shows several Weber waves for different values of the parabolic momentum. 



Note that the Weber waves exhibit a shape-preserving profile while moving on a parabolic 
trajectory whose shape is controlled by the parabolic momentum. The focus of the parabola 
increases as the parabolic momentum does. For a < the Weber waves will propagate to 
the left because W(x,z; —a) = W(— x, z; a). The bending angle 4>b (respect to the z-axis) 
of a Weber beam is given by <pb = (tt + 6)/2, where 9 is the polar coordinate angle, i.e., 
x = t cos 8 , z = rsin#. We found that the first two lobes of the "half a Bessel" and Weber 
waves at z = have the same shape if the parameters of the waves satisfy the relation 
m = 2a. In the top row of Fig. rtTj) we show several "half a Bessel" waves with m = 2a; 
although we can see that the "half a Bessel" bend faster, qualitatively, it is possible to see 
that Weber waves stay true to their propagation invariant nature longer. If we compare the 
Weber waves with "half a Bessel" waves that remain nondiffr active longer, for example the 
ones satisfying the relation m = 3a depicted in the bottom row of Fig. Q, one can see that 
the bending trajectories are equivalent. 

We found that for a > 10 our Weber waves can be approximated (with an error less than 
0.5% relative to the maximum) by 
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where x = X\fkj2a and 
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xV x 2 + 1 + log yx + V x 2 + 1 

This form of the Weber waves unravels their physical characteristics. Notice that now our 
waves separate in a product of functions of the r] and £ coordinates. The £ part describes 
the profile characteristics while the rj part describes the propagation characteristics. The £ 
part is a modulation of the Airy function, the argument of the Airy function makes this part 
oscillate asymptotically as sin (/c£ 2 /2 + 7r/4), i.e., the maximum frequency of oscillation of 
the wave will be k as expected because the Weber waves do not contain evanescent waves. 
While the modulation make the envelope of the function to decay as £ -1 ^ 2 , i.e., as an Airy 
beam. 
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FIG. 2. (Color online) The color map and parabolic lines shows the decay and propagation trajec- 
tory, respectively, of the main lobe of the Weber waves for different parabolic momentum a. The 
coordinates x and z are in units of k^ 1 . 

We found that the maximum lobe of the Weber beam propagates in the parabola with 

V^max « 2^a~ + 0.9838a" 4/25 , (9) 

and by Eq. the value of this maximum will decay upon propagation by 5 = (krf /2a + 1) ^ 4 . 
For Weber waves with a > 10 this decay can be simplified within an error of 5% by 
5 = a/cos 06 ; from this we see that our waves will bend to 75° with still 50% of its peak am- 
plitude. Figure [2] shows the decay and propagation trajectory of the main lobe of the Weber 
waves for different values of parabolic momentum, i.e., it characterizes the nondiffractive 
accelerating behavior of our waves. 

It is possible to double the angle of bending by propagating the Weber waves from z < 0. 
In this case the wave has a bending angle opposite to the direction of bending and the Weber 
wave will depict full parabolas as show in Fig. [3]^a). These Weber waves can be used for 
micromachining of micron-size curved structures [TO]. 

Because the Weber waves are separable solutions of the wave equation, they are also 
eigenfunctions of a second order operator that we will call the parabolic momentum operator 
|17j . This operator commutes with the Helmholtz equation; therefore, the eigenvalue of the 
parabolic momentum operator is a conserved dynamical constant. The parabolic momentum 
operator is given by 



(Ld z + d z L) W 
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where L = xd z — zd x = (r/d^ — (,d v ) /2 is the angular momentum operator. 
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FIG. 3. (Color online) (a) Weber wave a = 55; (b-e) Paraxial Weber waves with d = 3,7. (a) 
depicts a 440&; -1 x 220&; -1 rectangle and (b-e) all depict a 220&;~ 1 x 220&;~ 1 square. 

To explain the physics of the parabolic momentum let us calculate the local value of the 
linear momentum density of the Weber waves W = \W\ e l/3 , 

p = iuEo \W\ 2 V/3, 

= i ue \W\ 2 k^—J-rf- (11) 

Vv + r 

This means that the linear momentum flows purely in parabolas. Therefore, the trans- 
verse structure of the Weber waves allows trapping micrometer particles by gradient force, 
and the linear momentum density will move these particles in parabolic trajectories. This 
opens the possibility to apply schemes and uses equivalent to those found for orbital angular 
momentum [12] to the parabolic momentum. 

We will now compare the known features of the Airy beam and "half a Bessel" waves 
to our new Weber waves. First, the Weber waves are also self-healing, for example, in 
Fig. gja) the main lobe of a Weber wave is initially blocked and during propagation the 
second lobe gets more power and its trajectory bends to replace the first lobe, and similarly 
each subsequent lobe will bend to replace the lobe on its right. 

As in the case of the "half a Bessel" waves, by superposing Weber waves with different 
values of parabolic momentum, it is possible to create breather waves that will propagate in 
a parabolic trajectory but its amplitude will oscillate with periodicity that depends in the 
difference between the values of the parabolic momentum cf. Fig. Qb). Similar to the Airy 
beams and the "half a Bessel" waves our Weber waves are not square integrable, i.e., they 



carry infinite power. However, as it is common with accelerating beams, if we introduce an 
aperture or an exponential apodization, the Weber waves will hold their properties but over a 
smaller distance; the wider the aperture the larger the distance. Also, while the acceleration 
of the Airy beams can be controlled by an arbitrary constant that scales the linearly 
and the z— axis quadratically; the acceleration of our Weber beams is controlled by a the 
parabolic momentum, as we explained above. 

It is important to notice that although Airy beams are solutions to the paraxial equation 
they do not satisfy the paraxial approximation because its spectrum is not confined to small 
angles; this is equivalent to the fact that the Airy beams have an oscillating tail of constantly 
increasing frequency. Any solution of the wave equation only satisfies the paraxial equation 
in the paraxial regime; therefore, it is not possible to recover the Airy beams by varying the 
parameters of the "half a Bessel" wave or the Weber waves. Furthermore, neither the "half 
a Bessel" nor the Weber waves, in their current form, have a parameter that controls its 
paraxiality. In order to control the paraxiality of the waves we can use the idea of imaginary 
displacement [TBI [19] . If one considers a spatial shift by the imaginary distance id in the z di- 
rection, then the spectral function A{6) becomes A(6) exp(<isin#). Therefore, as d increases 
the spectrum becomes increasingly confined around the propagation direction and the field 
becomes increasingly beam-like. The amplitude of the spectral function of the Weber waves 
has a singularity at 9 = 0,7r, we will call this a soft singularity because as the amplitude 
diverge the phase also diverge; therefore, in any experimental or numerical application this 
divergence will average to zero. Nevertheless, it is possible to cancel this divergence if we 
subtract the original spectral function from the shifted one, i.e., A(8) (exp(kqsm8) — 1). In 
this case, the intensity of the spectral function at 9 = 0, it relative to the maximum inten- 
sity is d 2 / exp(2<i), consequently in any practical application we can consider these highly 
nonparaxial components zero if d > 3. Figure |3](b-e) depicts Weber waves constructed with 
this spectral function for d = 3, 7. 

Vector field solutions can be readily constructed by using our Weber waves as Hertz 
potentials to Maxwell equations, e.g., defining E TE = W(r],^)y and H™ = W(rj,^)y as in 
[TT] . Furthermore, arbitrary polarization vectorial Weber waves are given by Eq. (5) of 

To summarize, we have found nonparaxial accelerating waves that propagate along a 
parabolic trajectory while preserving its shape to a good approximation. These waves have a 
well-defined conserved dynamical constant that characterize them, the parabolic momentum. 
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FIG. 4. (Color online) (a) The self-healing property of a Weber wave (a = 40) with the first lobe 
initially cut out. The dashed black curved describes the original trajectory of the first lobe, (b) 
A periodic Weber wave created by supperposing two Weber waves with a = 40 and a = 54 with 
equal coefficientes. All figures depict a 220k- 1 x 220/c" 1 square. 

For moderate to large values of the parabolic momenta the Weber waves can be described by 
a modulated Airy function, therefore these waves by functional form are the closest related 
to the paraxial Airy beams. This work shows that there is another nonparaxial accelerating 
solution to the wave equation and that the "half a Bessel" waves are not an isolated case. 
This opens the possibility to explore new accelerating solutions of the Maxwell equations. 
Finally, because the Weber waves are exact time-harmonic solution of the wave equation, 
they have implications to many linear wave systems in nature, ranging from sound, elastic 
and surface waves to many kinds of classical waves. 
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